Reexamining the neutron skin thickness within a density dependent hadronic model 
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In the present work we calculate the 208 Pb neutron skin thickness, binding energy, surface energy 
and density profiles within the Dirac solution of a density dependent hadronic model. The same 
calculation is performed with the NL3 parametrization of the non-linear Walecka model. The asym- 
metry of a polarized electron scattered from a heavy target is also obtained within a partial wave 
expansion calculation. The theoretical results are then ready to be compared with the experimental 
results expected to be available very soon at the Jefferson Laboratory. For completeness, other 
nuclei as 40 Ca, 48 Ca, 66 Ni and 90 Zr are also investigated. 
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I. INTRODUCTION 

Relativistic models are a very useful tool in the de- 
scription of a wide variety of applications in nuclear mat- 
ter, finite nuclei and nuclear astrophysics. They can be 
tested and, hopefully constrained, according to exper- 
imental and astrophysical observation results. Unfor- 
tunately, most theoretical results are model dependent 
and so far it is unclear whether some of them should 
be discarded. Many variations of the well known quan- 
tum hadrodynamic model [l[ have been developed and 
used along the last decades. Some of them rely on den- 
sity dependent couplings between the baryons and the 
mesons 0, 0, 0, H| while others use constant couplings 
0) 0, HI- Still another possibility of including density 
dependence on the Lagrangian density is through the 
coupling of the mediator mesons among themselves 0. 
The strong model dependence of the results come from 
the simple fact that relativistic model couplings are ad- 
justed in order to fit expected nuclei properties such as 
binding energy, saturation density, compressibility and 
energy symmetry at saturation density only. Once the 
same relativistic models are extrapolated to higher den- 
sities as in stellar matter or higher temperatures as in 
heavy-ion collisions or even to lower densities as in the 
nuclear matter liquid gas phase transitions, they can and 
indeed provide different information. Hence, experimen- 
tal constraints obtained either from polarized electron 
scattered from a heavy target, from heavy- ion collisions 
at different energies or from astronomical observations 
are very important in order that adequate models are 
chosen and inadequate ones are ruled out. 

In the present work we focus our attention to the cal- 
culation of the difference between the neutron and the 
proton radii known as the neutron skin thickness. Based 
on the argument that experimental results should be used 
to constrain relativistic models, it is very important that 
an accurate experimental measurement of the neutron 
skin thickness is achieved. This depends on a precise 
measurement of both the charge and the neutron radius. 



The charge radius is already known within a precision of 
one percent for most stable nuclei, using the well-known 
single-arm and non-polarized elastic electron scattering 
technique as well as the spectroscopy of muonic atoms 
[l(J. For the neutron radius, our present knowledge has 
an uncertainty of about 0.2 fm However, using po- 
larized electron beams it is possible to obtain the neutron 
distribution, as first discussed in [l2j and, as a conse- 
quence, to obtain the desired neutron radius. In fact, the 
Parity Radius Experiment (PREX) at the Jefferson Lab- 
oratory [l3| is currently running to measure the 208 Pb 
neutron radius with an accuracy of less than 0.05 fm, 
using polarized electron scattering. In the PREX exper- 
iment, the asymmetry is expected to be measured at a 
momentum transfer q sa 0.4 fm -1 . 

At this point it is worth mentioning that two other 
experimental methods have been used in order to mea- 
sure neutron skins and neutron halos [3, [TH, namely, 
the nuclear spectroscopy analysis of the antiproton an- 
nihilation residues one mass unit lighter than the target 
mass and the measurements of strong-interaction effects 
on anti-protonic X rays. Whenever possible, the results 
obtained with the above mentioned methods are also used 
for comparison in the results section. 

In a recent work [lq |. the neutron skin thickness and 
the asymmetry for polarized electron scattering off a 
hadronic target were investigated with the help of both 
density dependent and with constant coupling relativistic 
models. For these calculations two important compro- 
mises were made: for the calculation of the proton and 
neutron densities, a Thomas-Fermi approach was used 
and for the calculation of the asymmetry a Plane Wave 
Born Approximation for the electron [171 ] was enforced. 
The 208 Pb neutron skin thickness was then obtained with 
two different density dependent hadronic (DDH) model 
parametrizations, the TW [3J, where the density depen- 
dence is introduced explicitly through the couplings and 
the NLwp model where the density dependence ap- 
pears through the coupling of the vector and the isovec- 
tor mesons. One of the most used parametrizations of 
the non- linear Walecka model (NLWM), the NL3 was 
also used. The asymmetry, in the momentum transfer 
range of interest for the calculation of neutron skins, was 
shown to give very similar results for all models. On the 
other hand, the neutron skin thickness is smaller with 
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the DDH model than with the NL3. As the coupling 
strength between the ui and p mesons increases in the 
NLwp model, the neutron skin thickness moves from the 
original NL3 towards the DDH results. As the momen- 
tum transfer increases, the asymmetry becomes strongly 
model dependent. 

In the present work we revisit the same quantities, but 
with improvements in both approximations mentioned 
above in order to check if the model differences remain. 
As the target is a heavy nucleus ( 208 Pb), the results for 
the asymmetry are reconsidered with the use of the par- 
tial wave expansion method as briefly discussed in the 
Appendix. On the other hand, the proton and neutron 
densities are calculated from the solution of the Dirac 
equation. We also extend our calculations to obtain the 
neutron skin thickness of 40 Ca, 4a Ca, 66 Ni and 90 Zr. For 
the asymmetry, we also include applications for 48 Ca. As 
the NLwp model was shown to give results that interpo- 
late between the NL3 parametrization of the NLWM and 
the TW parametrization of the DDH model, we next re- 
strict ourselves to the NL3 and the DDH models, i.e., 
one with constant couplings and another one with den- 
sity dependent couplings. We also comment on the dif- 
ferences between the results obtained in the present work 
and within the Thomas- Fermi approximation. It is worth 
mentioning that most applications to neutron stars, to 
equations of state used to describe supernova simula- 
tions and the description of nucleation processes have 
been done within the Thomas-Fermi approximation. Un- 
derstanding its limitations and accuracy is indeed very 
important. Information about the skin thickness as ob- 
tained in finite nuclei can be directly related to astrophys- 
ical calculations and so the possible connections between 
Thomas-Fermi and Dirac density profiles can be useful. 

The paper is organized as follows: in section II we show 
the Lagrangian density of the DDH model and describe 
the formalism used; in section III we present and discuss 
the results; in section IV we draw our final conclusions. 



II. FORMALISM 



where <f>, b M and are the scalar-isoscalar, vector- 
isoscalar and vector-isovector meson fields and the pho- 
ton field respectively, fi^„ = d^V v — d u V^ , — 
d^h v - <9„b M - r^b^ xb„), F^ v = d^A u - dyA^ and 
T3 = ±1 for protons and neutrons respectively. The pa- 
rameters of the model are: the nucleon mass M = 939 
MeV, the masses of the mesons m s , m v , m p , the elec- 
tromagnetic coupling constant e = y/ Aw/ 137 and the 
density dependent couplings T s , r„ and r p , which are 
adjusted in order to reproduce some of the nuclear mat- 
ter bulk properties, using the following parametrization: 



with 



and 



Ti (p) = Ti (p sat )hi (x) , X = p/p sa t , 



1 + hix + di) 2 
h i (x) = a i — — , , so , i=s,v 



1 + c^x + dif 



h p (x) 



exp[- 



1)], 



(2) 



(3) 



(4) 



with the values of the parameters m^, Ti(p sat ), a^, bi, 
Ci and di, i = s,v,p given in [3J]. This model does not 
include self-interaction terms for the meson fields as in 
NL3. We consider two parametrizations of the above 
mentioned DDH model, the original one that we next re- 
fer to as TW and a more recent one known as DDME1 
[HI , obtained from a fitting that includes known experi- 
mental values of the 20S Pb neutron skin. 

Once the Lagrangian density is chosen, the Euler- 
Lagrange equations are used to calculate the equations 
of motion. The meson field equations of motion are eas- 
ily found in the literature and we refrain from writing 
them. An interested reader can obtain the equations in 
0, [H, [lH, H3| , among other papers in the literature. 

The Dirac equation for the nucleon field reads: 



[y - E„) - (M - E,)]¥ = , 



(5) 



where the scalar and vector self-energies are given respec- 
tively by 



We describe the main quantities of the DDH model, 
which has density dependent coupling parameters in the 
following. The Lagrangian density reads: 



with 



E s = T^, E M - E£» + E* 



(6) 



= + • b M + e^A, , (7) 



+-(d^d^ - m 2 J 2 ) - -Jl^fF" 



\m 2 v V^ - ]b, v ■ B"" + \m 2 p h, ■ b^ - -F, V F^ (1) 



dp 



The term E? is known as the rearrangement term. As 
done in [l], , we consider that the nucleus is approxi- 
mately given by a Slater determinant and only the low- 
est lying positive energy states are occupied. The usual 
ansatz for the nucleon spinor is 



4>(j,t) = V>(r) exp(— iEt), 



(8) 
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and the one-particle states are obtained from the solution 
of the Dirac equation. Only spherically symmetric nuclei 
are considered and the usual notation for the expected 
values of the meson fields are used. The self-energies 
become: 

S = r„Vb(r) + ^r p r 3 b (r) + c (1 ± ^ A (r) + Eg , (9) 

S = [j^P(r)Vo + 5 ^P3(r)6o - ^P.(r)*J , 

(10) 

where P3 — p P ~ p n and p s and p are the usual scalar and 

baryonic densities. Next the ansatz given in eq.© and 

( 9 K (r)y 3m \ 
the spin-isospin wavcfunctions ib — [ . \*4 m <8> £ 

v *Mr)y- K J 

are substituted into the Dirac equation. y±^ and £ are 
the spinorial spherical harmonics and the isospin wave- 
functions respectively. After some straightforward alge- 
braic manipulations, a new system of coupled equations 
is obtained: 

(M*(r) + V(r)) g K {r) - ( £ - —) /„(r) = Eg K (r) 



It is also useful to define the nuclear charge distribution 
as: 

= [(Q 2 _ ffl )7r]3/a / d'r'expi- ^l^ ^'), 

(14) 

where B 2 = b 2 /A, a p — 0.653 fm and b is the oscillator 
length used to define our basis set. The above defini- 
tion includes the proton finite size, for which we have 
assumed a gaussian shape to describe its intrinsic charge 
distribution: 

Pint{r) = , 2 \ 3/2 exp (-(r/ap) 2 ), (15) 

and an approximate center-of-mass correction according 
to the prescription given in [2l| . Neutron finite size cor- 
rections are disregarded, i.e., its intrinsic form factor is 
taken as unit. 

For the calculation of the differential elastic electron 
cross-section and the asymmetry [A), details are 

given in the Appendix. 

III. RESULTS AND DISCUSSION 



d_ 

Or 



where 



and 



g K (r) - (M*(r) - V{r)) f K (r) = Ef K (r) , 

(11) 

M*{r) = M-r s (r)<f>{r), 



V(r) = T v (r)V (r)+-^T 3 b (r)+e 



A)(r)+E«(r). 



At this point, f K (r) and g K {r) are expanded in the har- 
monic oscillator basis of dimensions N and M respec- 
tively, as in The linear variational method is then 
used to solve the above equations and we are left with 
an eigenvalue problem solved by the diagonalization of a 
matrix of order N + M. The fields </>, Vb and bo are also 
expanded in a harmonic oscillator basis and solved iter- 
atively through their corresponding Klein-Gordon equa- 
tions. The Coulomb field Ao is obtained by the Green's 
function method. 

The proton and neutron mean-square radius are de- 
fined as 



J <Prpi[T) 
and the neutron skin thickness reads 

S = Rn — Rp- (13) 



We first plot the 208 Pb neutron and proton densities 
obtained for the DDH model within the Thomas-Fermi 
approximation from [l6l ] and the solution of the Dirac 
equation in Fig. [TJ As it is well known, the Thomas- 
Fermi approximation flattens the central energy densi- 
ties while the Dirac equation keeps the internal nuclear 
structure. The difference in the central density is then 
compensated near the surface, with consequences in the 
calculations of the mean square radii and neutron skin 
thickness. One can see that the DDME1 parametriza- 
tion gives almost the same density profiles as compared 
with the original TW parametrization. 

In table I we show the results of the present cal- 
culations for 208 Pb obtained with the NL3 and both 
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FIG. 1: Pb neutron and proton densities obtained with TW 
(Thomas-Fermi), DDMEl (Dirac) and TW (Dirac). Larger 
densities are for neutrons and lower for protons. 
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FIG. 2: Skin thickness as a function of the proton-neutron 
number asymmetry. Squares represent our values in the TW 
and full circles in the NL3. The lines represent the limits for 
the fitting given in [l4| ] . From lower to higher asymmetries we 
have plotted the results for 40 Ca, 90 Zr, 66 Ni, 48 Ca and 208 Pb. 

parametrizations of the DDH model. For the sake of 
completeness we also include the results obtained in [l6[ 
within the Thomas-Fermi approach. One can observe 
that the Thomas-Fermi results for the radii are system- 
atically larger than the Dirac solutions, but the neu- 
tron skins are smaller. The surface energies also give 
smaller contributions within the TF approach. Both 
parametrizations of the DDH model show the same neu- 
tron skins albeit the small differences in the proton and 
neutron radii. The binding energy and the surface en- 
ergy are also slightly different. Finally, the charge radius 
(-Rc)j calculated with the charge distribution given by 
equation (fT4f is presented and compared with the exper- 
imental value. All results are very similar. 

In table II we show the results for 40 Ca, 48 Ca, 6f Wi 
and 90 Zr. Most of the conclusions drawn from Table U 
remain valid, i.e., the neutron skin thickness are the same 
for both parametrizations of the DDH model, TW gives 
larger surface energies than DDME1 and NL3 and the 
proton, neutron and charge radii vary from one model to 
the other. 

According to [l4| . experimental data show a linear 
relation between the neutron skin thickness and the 
proton- neutron asymmetry of the considered nuclei, i.e., 
S = (N — Z)/A. Assuming that this dependence really 
exists, the authors of [3] fit it as 

6 = (-0.04 ± -0.03) + (1.01 ± 0.15)5. (16) 

In Fig. [5]we have plotted the fitted dependence lines con- 
sidering the most extreme values for 0. If this dependence 
has really to be satisfied, the theoretical points should lie 
within both curves. The TW and DDME1 points co- 
incide and they are shown to be inside the appropriate 
range while the NL3 results for very asymmetric nuclei 
are outside the upper boundary. 

Now let's focus on the electron scattering cross sec- 
tions and the asymmetry. In Fig. [3] the elastic cross 
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FIG. 3: Elastic electron scattering cross section in the NL3 
(dashed line), TW (full line) and DDME1 (dotted line) ap- 
proaches for a) 208 Pb and b) 48 Ca. E is the electron incident 
energy taken as 502 MeV for 208 Pb and as 757.5 MeV for 
48 Ca. Experimental data are from [221 ] and [23I ] respectively. 
The dot-dashed line in a) corresponds to the Thomas-Fermi 
approximation and TW parametrization. 

section for NL3 and DDH models are compared with the 
experimental data for 208 Pb and 48 Ca. Finite size pro- 
ton as well as center-of-mass effects were also included in 
the cross-section calculation. In both cases the theoreti- 
cal and experimental results are in very good agreement 
and both models show very small quantitative differences 
along all the momentum transfer region covered by the 
data. While TW and DDME1 seem to give a slightly 
better desciption for 208 Pb, NL3 results are closer to 
the cross section for 48 Ca at momentum transfers be- 
tween 3.0 and 3.5 fin -1 . For the sake of completeness 
we have added the results obtained with the Thomas- 
Fermi approximation for the TW parametrization of the 
DDH model for 208 Pb. One can see that the results for 
the cross section start to deviate from the experimen- 
tal points around 1.5 fm _1 with larger discrepancies at 
larger momentum transfers. 

As for the asymmetry, our results are displayed in Fig. 
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IV. CONCLUSIONS 

In the present work we have recalculated the neutron 
skin thickness and the asymmetry that will be measured 
at the PREX experiment with two kinds of relativis- 
tic models previously used within a Thomas-Fermi plus 
PWBA approximations. We have here considered the full 
solution of the Dirac equation and used an exact calcu- 
lation for the scattered electron wavefunction. Density 
dependent and constant coupling relativistic models pro- 
vide different results and so the model dependence of the 
electron scattering asymmetry is confirmed by the more 
exact calculation, although it is still very hard to be ex- 
tracted from small momentum transfer data. 

If a linear relation between the neutron skin thick- 
ness and the proton-neutron asymmetry of the consid- 
ered nuclei is really to be satisfied, as suggested in [Til ], 
the TW and DDME1 model parametrizations provide re- 
sults within the appropriate range while the NL3 results 
for very asymmetric nuclei are outside the upper limit 
imposed by the present data. 

Once a more precise measurement of the neutron skin 
is published, it will certainly provide a better constraint 
to different models and our results should then be revis- 
ited for a proper comparison. 



FIG. 4: Asymmetry obtained for a) 208 Pb and b) 48 Co in the 
NL3 (dashed line) and TW (full line). The dotted line is the 
PWBA result and the dash-dotted one is the three parameter 
Fermi result (FermiSp) as explained in the text. The incident 
electron energy was chosen as E = 800MeV\ 



@]for 208 Pb. The curve labeled PWBA is simply the re- 
sult obtained from equation (|23|) in the Appendix, for 
the case where Z p n = N p p . This can be compared 
with the curve labeled FermiSp, where the same condi- 
tion is imposed but the cross-section calculation consid- 
ered the partial wave expansion method and a three pa- 
rameter Fermi distribution to obtain the Coulomb and 
weak potentials as defined in the Appendix. The dif- 
ference between those curves clearly states the necessity 
to incorporate the electron wave distortion effects in the 
calculation of the asymmetry. The other two curves give 
us the structure effects coming from the NL3 and DDH 
models respectively. Similar results are shown in Fig. [4] 
for 48 Cet. One should notice that while the proton dis- 
tribution in the NL3 and DDH models follow the same 
trend even for large momentum transfer, as can be seen 
from Fig. [3j the asymmetry starts to present important 
qualitative differences at q > 2/to~ 1 , indicating that the 
neutron distribution seems to be more sensitive to the 
model used. 



V. APPENDIX - CROSS SECTION FOR 
POLARIZED ELECTRONS 

The electron hadron Lagrangian density is 

C = $ bv (i0" + eA" + - m e ] V + C H , (17) 

where Ch is the hadronic Lagrangian, is the electro- 
magnetic field generated by the hadronic target, 



NC 



(18) 



is the weak held and G the Fermi constant. In the elastic 
scattering of electrons from an even-even spherical nu- 
cleus whose elementary particles are just the protons and 
neutrons, only the static time component of the currents 
contribute and the electron obeys the Dirac equation: 



[a.~p + 7 m e + V(r) + j 5 V A (r)]ip = Eip, (19) 

with V(r) = -eA {r) and V A {r) = ^J 7VC (r). From 

the standard model we know that the weak neutral (NC) 
current is 281: 



J°' NC (r) = x P Pp{r) + XnPn(r). 



(20) 



The above constants are Xp — 0-04 and Xn = —0.5. 
Hence, we conclude that the weak potential depends 
strongly on the neutron distribution in nuclei, while the 
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cm. field Aq depends (not considering the neutron elec- 
tric form factor) exclusively on the proton distribution. 
For high energy electrons and scattering angles not to 
close to 180° we may use the approximation ^ ~ 0, for 
which the Dirac equation can be rewritten in the form: 



polarized electrons with positive(+) and negative (— ) he- 
licities. The solution of equation (f2"Tj) in first order per- 
turbation theory lead us to the Plane Wave Born Ap- 
proximation [12| result: 



. + V{r) ± VA(r)]V± =£V±. 



(21) 



with tp± = ^(1 ± 75)^ and the ± signs represent the two 
possible electron initial polarization (hclicity) states. We 
have solved the above equation for the electron exactly, 
using the well-known partial wave phase-shift expansion 
method, as explained, for instance in [29]. Note that for 
each electron helicity state, a different set of phase-shifts 
must be determined, whose differences come from the 
contribution of the weak potential. We now define the 
asymmetry through the expression: 



A 



Gq 2 



\Xp + X; 



Pn(g) i 
'Pp(l) 



(23) 



where q is the momentum transfer and p n ( p ){q) being 
the neutron (proton) distribution in q space. Even for 
medium mass nuclei, important differences between the 
PWBA and the exact result should be found in the asym- 
metry [T3|. In this paper, all the cross sections were ob- 
tained through the full solution of the Dirac equation for 
the electron. 



da+ I d£l — da~- 1 dil , „ 
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TABLE I: 208 Pb properties 



model 


approximation R n R p R c 


e 


B/A 


a 




(fm) (fm) (fm) 


(fm) 


MeV Mev/fm 2 


NL3 


TF 5.79 5.57 


0.22 


-7.79 


0.96 


NL3 


Dirac 5.74 5.46 5.51 


0.28 


-7.91 


1.13 


TW 


TF 5.68 5.52 


0.16 


-7.46 


1.10 


TW 


Dirac 5.61 5.42 5.48 


0.20 


-7.78 


1.30 


DDME1 


Dirac 5.66 5.46 5.51 


0.20 


-7.91 


1.18 


exp. [24] 


5.50 








exp. [25] 






-7.87 




exp. [26] 




0.12 ±0.07 






exp. [27] 




0.20 ± 0.04 






exp. [15] 




0.16 ±0.02 ±0.04 






TABLE II: Finite nuclei properties 



model nuclei R„ R p R c 9 B/A a 
(fm) (fm) (fm) (fm) MeV Mev/fm 2 



NL3 40 Ca 3.32 3.37 3.43 -0.05 -8.62 1.48 

TW 40 Ca 3.28 3.33 3.39 -0.05 -8.36 1.60 

DDME1 40 Ca 3.32 3.37 3.43 -0.05 -8.62 1.45 

exp. [18] 40 Ca 3.48 -8.55 

NL3 48 Ca 3.60 3.37 3.44 0.23 -8.72 1.54 

TW 48 Ca 3.54 3.35 3.42 0.19 -8.49 1.70 

DDME1 48 Ca 3.58 3.39 3.46 0.19 -8.66 1.53 

exp.[30] 48 Ca -8.67 

exp. [18] 48 Ca SA8 

NL3 90 Zr 4.30 4.19 4.25 0.11 -8.86 1.37 

TW 90 Zr 4.24 4.15 4.22 0.08 -8.55 1.52 

DDME1 90 Zr 4.28 4.19 4.25 0.08 -8.73 1.38 

exp. [30] 90 Zr -8.71 

NL3 66 Ni 3.96 3.76 3.82 0.20 -8.74 1.47 

TW 66 Ni 3.89 3.74 3.81 0.15 -8.56 1.63 

DDME1 66 Ni 3.93 3.77 3.84 0.16 -8.72 1.49 

exp. [30] 66 Ni -8.74 



